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In generic holographic QCD, we find that baryons are bound to form a nucleus, and that its 
radius obeys the empirically-known mass number {A) dependence r oc A^^'^ for large A. Our 
result is robust, since we use only a generic property of D-brane actions in string theory. We also 
show that nucleons are bound completely in a finite volume. Furthermore, employing a concrete 
holographic model (derived by Hashimoto, lizuka, and Yi, describing a multi-baryon system in the 
Sakai-Sugimoto model), the nuclear radius is evaluated as 0(1) x A^^^ [fmli which is consistent with 
experiments. 



To describe atomic nuclei directly by strongly coupled 
quark dynamics, QCD, is a long-standing problem in nu- 
clear physics and particle physics. It was only recent that 
lattice QCD simulations reproduce qualitatively the nu- 
clear forces. Recent progress in solving strongly coupled 
gauge theories with a new mathematical tool of super- 
string theory, the AdS/CFT correspondence, has been 
proven to be truly powerful in application to QCD (called 
holographic QCD). 

In this letter, we show by quite a generic argument 
in superstring theory and the AdS/CFT correspondence 
that non-supersymmetric QCD-like theories in the large 
Nc limit host nuclei, multi-baryon bound states. Fur- 
thermore, we can show also that the resultant nuclei have 
the important nuclear property in the real world: Finite- 
ness of the nuclear size, and its mass-number dependence. 
That is, the radius of the holographically realized nuclei 
is shown to be proportional to A^^^ where A is the mass 
number (the baryon number) of the nucleus. 

In deriving these, we do not rely on any specific model 
of holographic QCD. What we use is only the following 
two known facts: (i) Baryons are D-branes in any gravity 
dual of QCD-like gauge theories [21 E], (ii) D-brane effec- 
tive actions are a dimensionally reduced Yang-Mills (YM) 
theory [3]. From these two, the formation of the nuclei 
and the mass-number dependence of the nuclear size fol- 
low. Therefore our finding is quite robust and universal 
for any holographic description of non-supersymmetric 
QCD-like gauge theories, at the large Nc and at the 
strong coupling. 

Our derivation is divided into two steps. 
1. The system with a large number of baryons in 
generic holographic QCD is described effectively by 
a simple bosonic matrix quantum mechanics. It is 
a pure YM action dimensionally reduced to 1 di- 
mension 1261. 



dt iiA 



{DtX'f + Ux',X'] 



(1) 



where / = 1, • • • , D. The eigenvalues of the Ay. A 
matrix X* (i = 1, 2, 3) are location of the A baryons 
in our space, and X* (i = 4, • • • , _D) is for holo- 
graphic directions. 



2. The system allows a non-perturbative vacuum at 
which the eigenvalues of form a ball-like distri- 
bution, which is nothing but a nucleus. The size 
shows the mass-number dependence A^/"^. 

In the following, we shall show 1 and 2 in turn. Fi- 
nally we present the explicit form of the nuclear den- 
sity distribution (10). Together with the explicit matrix 



model [T] where input parameters are only the p meson 
mass and the NNtt coupling, we obtain the nuclear ra- 
dius R ~ 0{1) X A^/"^ [fm], which is consistent with the 
standard experimental observation. 

Multi-baryons described by matrices. — First, let 
us show that the system of A baryons in generic holo- 
graphic QCD is dictated by the matrix quantum mechan- 
ics ([T]). This can be provided by known facts (i) and (ii) 
which we explain below. 

(i) Baryons are D-branes in any gravity dual of QCD- 
like gauge theories, basically because baryons in large Nc 
gauge theories are heavy as their mass diverges as 0{Nc) 
while D-branes are solitons of string theory and thus 
heavy. In fact, from the point of view of baryon charges, 
it is known that, in the gravity dual side of strongly cou- 
pled gauge theories, D-branes wrapping compact cycles 
can be identified as baryons jU [3] ■ This is natural since 
in large Nc QCD baryons appear as solitons of meson 
effective field theory, as in the famous Skyrme model. 

To show the universality of the statement, let us illus- 
trate this fact with a specific and concrete dual gravity 
background provided by Witten [6| . The spacetime ends 
smoothly at an IR end capped with non-contractable S"*. 
A D-brane corresponding to a baryon is the one wrapping 
it (note that it is orthogonal to our spacetime so the D- 
brane looks as a point particle in our space) . As the is 
supported by Nc units of a Ramond-Ramond (RR) flux 
in the supergravity solution of the geometry, the wrapped 
D-brane, called a baryon vertex, has to attach the ends 
of Nc fundamental strings, thus attain a single baryon 
charge (= Nc quark charge) of the gauge theory [51. This 
situation is universal in holographic QCD, by the follow- 
ing two reasons. First, it is expected that, in any grav- 
ity dual of confining gauge theories, the geometry pos- 
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FIG. 1: Gravity dual of QCD-like confining gauge theories with 
A baryons. The horizontal directions are our 3-dimensional space, 
and the vertical direction is the holographic direction. This 10- 
dimensional spacetime is curved, and the bottom is capped by a 
smooth compact manifold, which is the IR end. The baryon vertices 
{A red balls) sit on the Nj flavor D-branes. The degrees of freedom 
of the baryons are open strings attached to the baryon vertices. 

sesses a smooth IR end accompanying a compact man- 
ifold, known as the name of "confining geometry" |27j. 
This end is necessary to provide a finite QCD string ten- 
sion. Second, the compact manifold always carries the 
Nc RR flux because the geometry is a gravity dual of 
SU{Nc) gauge theories realized by TV^ D-branes. 

Any nucleus is labeled by spins and isospins, so we 
need to assign flavor (isospin) charges on the baryon in 
the AdS/CFT. Here we review the standard holographic 
description of the flavor symmetry. The quark sector 
of QCD-like gauge theories is realized in the AdS/CFT 
correspondence by introducing flavor D-branes [5]. Typ- 
ically we need Nf D-branes which extend also along our 
spacetime directions, where Nf is the number of quark 
flavors. The iVj, fundamental string emanating from the 
baryon vertex end on these flavor branes. The excita- 
tion of these strings represent the flavor charges of the 
baryons. This is nicely realized in Sakai-Sugimoto model 
of holographic QCD [9]. Quantum excitations of the 
strings give rise to the spins and isospins [T] [3S|. In 
Fig. [T] we show a typical brane configuration of gravity 
dual of confining gauge theories with A baryons. 

(ii) Effective theory on the A baryon vertices for large 
A is given by the matrix quantum mechanics Q. Funda- 
mental degrees of freedom on any collection of D-branes 
are strings connecting the D-branes. In our case, we have 
A baryon vertices and Nf flavor D-branes, thus we have 
Ax A matrix from strings among the A D-branes, in 
addition to A x Nf matrix w coming from strings con- 
necting the baryon vertex and the flavor D-brane. These 
strings are shown in Fig. [l] The effective theory is a 
U{A) gauge theory, where is in the adjoint represen- 
tation while w is in the fundamental representation. The 
index / runs for our spatial directions ? = 1, 2, 3, and the 
holographic directions i = 4, • • • , D. 

Our interest is a large A, that is, a heavy nucleus, 
as A is the mass number. At the leading order in the 
large A expansion, the fundamental representation field 



w drops off. It is known in string theory that the X^ part 
has a universal effective action which is a dimensionally 
reduced YM action, ([T]). 

We ignore the gauge field components along the com- 
pact manifold {S'^ in the previous example), since those 
directions are irrelevant to QCD The time com- 

ponent gauge field At still remains, but it is not dy- 
namical and just ensures the local gauge invariance [30] . 
Fermionic superpartners are heavy due to the supersym- 
metry breaking and assumed not effective. 

An explicit action can be derived once we fix the 
species of the branes in holographic QCD. The nuclear 
matrix model of [1] [31] was derived in the Witten's ge- 
ometry with the fiavor D8-branes of the Sakai-Sugimoto 
model. 



S^^^Jdttr^^^{DtXr + ^^[X',X-'f 



-^M^kHX ) + ) ) j + sub-leading in l/A. (2) 

Here / runs from 1 to = 5, A = -^cffqco ^^^^ QCD 
'tHooft coupling, and Mkk is a dynamical scale (roughly 
corresponding to the QCD scale). With a rescaling = 
X^ {XNcAIkk/'^^t^Y^^ , the leading order is written in a 
canonical expression. 



S ^ dt tr 



i(AF^)^-^(r^) 



i\2 



+ ^[y^W]^),(3) 



where the matrix model coupling squared is = 
2^AM|^j^/(3^7rA^c)- In comparison to ([T]) here there is a 
mass term with my = (2/3)M^j^, but it is qualitatively 
irrelevant in our analysis, see Footnote [32]. The mass 
term is only along the holographic direction X^, which 
means that the baryon vertex is stable at the IR end of 
the geometry X"^ — 0. 

Formation of nucleus and nuclear size. — We shall 
show that indeed the eigenvalues of the matrix X^ are 
bound to each other, which directly means a formation 
of a nucleus, since the eigenvalues are the location of the 
A baryons. 

As the action ([T]) is a dimensional reduction of a YM 
theory, we can apply the generic argument by Liischer 
[17] on spectra of YM theory on torii. Liischer proved 
a theorem that the ground state of the theory ([I]) has 
eigenvalues bound to a finite region, and is invariant un- 
der the U{A) rotation and the spatial rotation. Even 
though there exists a flat direction of the potential at 
which all X^ are diagonal, the eigenvalues are bound be- 
cause the flat direction is narrow for large values of X^ 
and the quantum dynamics suppresses percolation. So, 
this theorem ensures a formation of a spherical nucleus 
in generic holographic QCD. 

Next, we show the nuclear radius oc It follows 

simply from a dimensional analysis of the action ([ij . The 



nuclear size is given by the distribution of the entries of 
(i = 1, 2, 3), so we can define the mean square radius 

as 



2 ^ 
^mcan — a. 



1=1 



(4) 



The expectation value is taken with the ground state of 
the matrix quantum mechanics. The normalization \/A 
is understood from the case of diagonal X with which 
individual location of the baryons makes sense. Using 
the properly normalized action as in (|3| , we can estimate 



^(tr^[rr])=coA; 



-1/3 



-o(iM), 



(5) 



at large A, where Xa = Ag^ is a 'tHooft coupling of 
the quantum mechanics and cq is an undetermined A- 
independent dimensionless constant. We have used a 
standard 'tHooft expansion [TH]. The Xa dependence 
was determined from the mass dimension: [Xa] = 3 and 
[y*] = —1/2. {Xa is the unique dimension-ful parame- 
ter in the action.) Since the rescaling from X to F is 
A-independent, we obtain fmean A^^^. Therefore 
we conclude that in generic holographic QCD a nucleus 
forms as a bound state of A baryons at large A, and ex- 
hibits the correct A dependence of the nuclear size. 

Nuclear density distribution. — A certain approxi- 
mation of the matrix quantum mechanics ([T]) enables us 
to compute more detailed information of the bound state: 
nuclear density distribution. Here we shall show that the 
nuclear density computed by the matrix quantum me- 
chanics ([T]) vanishes outside a certain radius, which is 
the finiteness of nuclei. 

To this end, we employ a RR density formula [19 in 
a D dimensional spacetime, at the leading order in 1/A, 
developed in the context of Matrix theory [20] in super- 
string theory: 



p{x) 



(2^) 



D 



d^fce-**-^ {tTA exp[ik-X]) 



(6) 



Since our baryon vertices are RR-charged D-branes, ^ 
is equivalent to the baryon charge distribution. We eval- 
uate the expectation value at zero temperature |32) . 

To evaluate ([6| with the model Q, we need a help of 
a large D expansion of matrix models pTH25] . where D 
counts the number of the matrices X^ (/ = 1, • • • ,D). 
We fix Xa = <?o^-D to be finite, and take the large D and 
large A limits. The large D limit is known to be a good 
approximation as discussed in [35] even for small D (> 2) 
qualitatively. According to [Hj, at the leading order of 
the expansion, we obtain a non-perturbative vacuum at 
zero temperature, which is characterized by 



-if (tr^[r^y^]),^„ = ^A-/^ 
1=1 



(7) 



p(r) 




ro r 

FIG. 2: The nuclear density distribu- 
tion p(r) obtained in the large D limit, 
ITo) . At the core of the nucleus r ^ 0, 
the distribution is homogeneous. For 
r > ro, the density is exactly zero, 
meaning the complete finiteness of nu- 
clei. A graphical image is drawn in the 
next figure. 



FIG. 3: The density 
distribution of the nu- 
cleus. Darker color 
means smaller density. 
Ignoring the surface 
defect, we find a homo- 
geneous density distri- 
bution. 



This indeed is consistent with Around this vacuum. 



1/3 
A ' 



as 



behaves as a free massive field with a mass A 
the interaction is suppressed by [33 . By using the 
propagator of the free massive scalar at zero temperature, 



{Ylbm/M)T=0 ~ -J^SadSbcS" {t ^ 0), 

2A^ 



(8) 



? can evaluate ^ as 

oo 

(tr^ exp(ifc • Y))t=o = ^ 



-A^ 



= A- 



ro\k 



n\{n + 1)1 \ 2X]/^ 
' Ji('^olfcl) 



(9) 



where tq = {2A/X]('^y^''^. Here we have used a fact that 
only ladder diagrams [23] contribute to ([6]). The distri- 
bution in our 3-dimensional space can be obtained by 
a Fourier transform followed by an integration over the 
holographic directions, 



D 



p{r) 



Y[dx- 

i=4 



1 



(27r)^ 
A 



= { T^'^-rl^Jrl - r2 

(r-o < r) 



d^k e-^^-''{XTA exp(ifc • Y))t=^ 
(r < ro) 



(10) 



where r is a radial coordinate in the 3-dimensional space 
spanned by {Y^, Y^,Y^}. We found that nuclear density 
is zero outside a cetrain radius, which is the finiteness of 
nuclei. 

Plotting this function shows the density distribution 
given in Fig. [2| and Fig. [3] We notice that at r = tq the 
density goes up, which is not really the case for realis- 
tic nuclei. However, since this surface part is at a sub- 
leading in the large A expansion, explicit computations 
of 1/A corrections may be necessary to see the details. 
Other possible origin may be the formula ^ itself [25j . 
We hope to come back to this issue in a near future. 
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We also notice that the isospin structure of the nu- 
cleus is invisible, since it is encoded in w, which is sub- 
dominant in the large A expansion. In contrast, the dis- 
tribution in X* directions indicates that the nucleus in- 
cludes components of excited baryon resonances. 

Evaluation of nuclear size. — Finally we shall study 
numerics. Using the explicit matrix quantum mechanics 
([2]) of jl], we can calculate Q through ^ as 



T=0 



35/2^2/3 ^1/3 
25/651/6 MKK{NcX^y/^' 



(11) 



A natural choice of the experimental inputs should be 
g-TrNN ^ 13.2 and uip ~ 776 [MeV] as these dictate nu- 
clear forces (the model [T] is in the chiral limit). From 
the work 13' we find that these experimental inputs are 
reproduced when A = 5.31 and Mkk = 949 [MeV]. This 
input leads us to a nuclear mean square radius 



V ' mi 



0.7 A^/^ [fml 



(12) 



As an order estimate, this is quite close to the observed 
value \/r^^^ — l.OA^^^ [fm] (which is equivalent to the 
well-known nuclear radius R ~ 1.2A^^^ [fm] as rf^^can = 
(3/5)i?^ for a homogeneous nuclear density.) 

One can alternatively think of tq in the density distri- 



bution ( 10 1 as the nuclear radius. After going back from 
Y to X, with the same experimental values as the inputs, 
we obtain tq ^ 0.8 x yl^/"^ [fm]. This again is consistent 
with experiments. 
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